A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

12(02. 01. 2;02.13.2;08.06.2;09.10.1;11.14.1;11. 10.1) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



Magnetically- driven jets from Keplerian accretion discs 



Jonathan Ferreira 

Landessternwarte, Konigstuhl, D-69117 Heidelberg, Germany 
jferreir@lsw.uni-heidelberg.de 

Received March 30; accepted July 2, 1996 



ON 
ON 



> 

in 
o 
i> 
o 

\o 

OS 
Oh; 

6 



Abstract. Non-relativistic, magnetically-driven jets are 
constructed by taking sclf-consistently into account the 
feedback on the underlying accretion disc. It is shown that 
such jets are mostly described by the ejection index ^ = 
d In Ma/d Inr, which is a local measure of the disc ejection 
efficiency. This parameter is found to lie in a very narrow 
range, due to constraints imposed by both the disc vertical 
equilibrium and the steady transfer of angular momentum. 

The investigation of global disc-jets solutions provided 
two important results. First, it shows that the disc vertical 
equilibrium imposes a minimum mass flux ejected. Thus, 
one cannot construct jet models with arbitrarily small 
mass loads. Second, their asymptotic behaviour critically 
depends on a fastness parameter uja — ^*rA/VAp,A, ratio 
of the field lines rotation velocity to the poloidal Alfven 
velocity at the Alfven surface. This parameter must be 
bigger than, but of the order of, unity. 

Self-similar jets from Keplerian discs, after widening 
up to a maximum radius whose value increases with uja, 
always recoUimate towards the jet axis, until the fast- 
magnetosonic critical point is reached. It is doubtful that 
such solutions could steadily cross this last point, the jet 
either ending there or rebouncing. Recollimation takes 
place because of the increasing effect of magnetic con- 
striction. This systematic behaviour is due to the large 
opening of the magnetic surfaces, leading to such an efh- 
cient acceleration that matter always reaches its maximum 
poloidal velocity. This "over-widening" stems from having 
the same ejection efficiency ^ in the whole jet. 

Realistic jets, fed with ejection indices varying from 
one magnetic surface to the other, would not undergo rec- 
ollimation, allowing either cylindrical or parabolic asymp- 
totic coUimation. The study of such jets requires full 2-D 
numerical simulations, with proper boundary conditions 
at the disc surface. 

Key words: Accretion, accretion discs - Magnetohydro- 
dynamics (MHD) - Stars: formation - ISM: jets and out- 
flows - Galaxies: nuclei - Galaxies: jets 



1. Introduction 

It is now widely accepted that accretion of matter onto a 
central object plays a major role in astrophysics: in active 
galactic nuclei (AGN), young stellar objects (YSOs), but 
also in evolved binary systems like cataclysmic variables 
or low-mass X-ray binaries. 

The standard theory of accretion discs (Shakura & 
Sunyaev |1973|, Novikov & Thorne |1973|, Lynden-Bell & 



Pringle 1974) assumes the presence of a turbulent viscos- 
ity that allows the matter to loose its angular momentum 
and mechanical energy. This energy is then radiated away 
at the disc surfaces, providing the observed luminosity 
(e.g., Bregman 1990t Bertout et al. 1988 ). However, such 



a theory is unable to explain the origin, acceleration and 
collimation of bipolar jets, that are emanating from radio 
loud AGN and quasars (Bridle & Perley |1984| ), all YSOs 
(Lada 1985 ) and some galactic objects like SS 433, Sco 
X- 1 (P ad man et al. 1991) and "microquasars" (Mirabel et 



al. 119921 ). 



In the vicinity of a black hole, the jet plasma has 
to come from a surrounding accretion disc. Lynden-Bell 
( 1978) suggested that in the framework of a thick accretion 
disc, jets could be produced in the inner funnel around the 
axis of the torus, and accelerated like in a De Laval noz- 
zle. Such a thick torus could be either radiation-supported 
(with a super-Eddington accretion rate, e.g. Abramowicz 
et al. 1980), or ion-supported (with a sub-critical accre- 



tion rate, Rees et al. 1982). Jets would then be accelerated 
by radiation pressure (Abramowicz & Piran 198C) or by 
the rotational energy of the central black hole (Blandford 
& Znajek 1977). This latter process, unique possibility 



within the second scenario, invokes a large scale magnetic 
field dragged in by the disc, braking the lastly rotating 
hole and transferring its energy to matter (see also Phin- 
ney 1983 ). However, it has been shown that such thick 
tori are viole ntly u nstable (Papaloi zou fc Pringle 1984 , 
Zurek & Benz 1986 , Begelman et al. |l987 ), casting there- 
fore strong doubts upon their viability to produce jets. 

The remaining scenario invokes a large scale magnetic 
field, anchored on a geometrically thin (Keplerian) ac- 
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cretion disc (Blandford & Payne |1984 hereafter BP82). 
These authors showed that this magnetic field could brake 
the disc and carry away its angular momentum. If jets 
carry a current, then the toroidal component of the mag- 
netic field would provide a tension that naturally confines 
the jet (previously recognized by Chan & Henriksen [19801 ). 

For protostars, the situation is more complex, since 
the observed jets could be either stellar winds (Canto 
1980| , Hartmann et al. 1982, Lago 1984) or disc winds 



(Pudritz & Norman [1983| Uchida & Shibata |l98|) . How- 
ever, both mass and momentum fluxes of the observed 
outflows are much higher than the ones provided by the 
protostar luminosity, hence forbidding both thermally and 
radiation pressure driven stellar winds (DeCampli 198l| , 
Konigl 1986). The possibility remains that stellar mag- 



netic fields play a major role in producing a wind (Mes- 
tel |1968[ Hartmann & Mac Gregor |1982|, Sakurai |198 



Tsinganos & Trussoni 1991, Sauty & Tsinganos |199 



But the strongest argument in favour of disc-driven jets 
is certainly the observed correlation betwe en signatures of 
accretion and ejection (Cabrit et al. 1990, Hartigan et al. 



19951) 



To summary, one can conclude that both observational 
and theoretical investigations tend to show that, in or- 
der to produce powerful self-coUimated jets, one has to 
rely on an accretion disc and a large scale magnetic field. 
Of course, the possibility that jets arise from the interac- 
tion between the central object magnetosphe re and the 
disc has st ill to be worked out (Camenzind |l990 Shu 



et al. 1994 ). Nevertheless, the study of such an interac- 



tion requires first the deep understanding of the interplay 
between accretion and ejection processes. Hereafter, we 
call Magnetized Accretion- Ejection Structures (MAES), 
objects where these two processes are interdependent. 
There have been a number of studies of magnetized jets 



in the past (BP82, Camenzind 1986, Lovelace et al. 1987, 
Heyvaerts fc No rman |1989|, C hiueh et al. |1991| , P ellcticr 
& Pudritz |1992t Li et al . [199^ , Appl & Camenzind |1993a, 



Rosso & Pelletier 1994, Contopoulos 1995, to cite only 



a few), but all these works were focused on jet dynamics. 
Thus, it was not yet proved that the underlying disc could 
indeed provide the required boundary conditions. 

Despite serious advances in the theory of magnetized 



accretion discs driving jets (Konigl 1989, Ferreira & Pel- 
letier [19934 |l993b| [1995| , Wardle & Konigl |1 993] , Li [1995] ), 
this question has not yet been fully addressed. Indeed, in 
both Wardle & Konigl and Li approaches, the disc so- 
lutions were obtained by not properly treating the disc 
vertical equilibrium and were directly matched to BP82's 
jet solutions (see Sect. 4.2.1). Moreover, by doing so, 
they were not able to specify the physical process lead- 
ing plasma to change its radial motion (accre tion) into a 
vertical one (ejection). Ferreira & Pelletier ( 1995 , here- 
after FP95) constructed disc solutions by taking into ac- 
count all the dynamical terms, thus being able to answer 
this question, as well as derive the physical conditions re- 



quired to magnetically launch jets. However, it remained 
to be proved that their solutions could indeed produce 
super- Alfvenic jets. 

The goal of this paper is therefore to construct global, 
non-relativistic solutions for magnetically-driven jets from 
Keplerian discs. Our treatment allows a smooth transition 
between the resistive disc and the ideal MHD jet. This pa- 
per is organized as follows. We start by briefly recalling 
the MHD equations of MAES and their main features. In 
particular, we expose some general results on the physical 
processes that govern the disc. Section 3 is devoted to the 
properties of magnetically-driven jets and the constraints 
arising from their interplay with the underlying disc. Self- 
similar, non-relativistic solutions that smoothly cross both 
slow-magnetosonic and Alfven critical points are displayed 
in Section 4. In the following section, we show the asymp- 
totic behaviour of self-similar jets and investigate analyt- 
ically the reason for their systematic behaviour. We con- 
clude by summarizing our results in Section 6. 

2. Keplerian accretion discs driving jets 

2.1. Magnetohydrodynamic equations for MAES 

In order to produce bipolar jets, it is natural to rely on 
a bipolar topology for the large scale magnetic field. The 
other alternative, a quadrupolar topology, could in prin- 
ciple be used too, but such a topology could only produce 
weak jets (if any, see Appendix ^). In both cases, this 
large scale magnetic field has two distinct possible ori- 
gins: advcction of interstellar magnetic field with accret- 
ing matt er or "in-situ" production th rough a disc dynamo 
(Pudritz 1981, Khanna & Camenzind 1994). Both scenarii 
raise unanswered questions. For example, on the degree 
of diffusion of the infalling matter, which strongly deter- 
mines the field strength in the inner regions (Mouschovias 
1991 , FP95). On the other hand, dynamo theory is still 
kinematic and one cannot easily infer from these studies 
what would be the final stage of the magnetic topology, 
if the matter feedback on the field is taken into account 
(Yoshizawa & Yokoi 1993). Most probably, a realistic sce- 



nario would have to take into account advection of ex- 
ternal magnetic field while its amplification by the local 
dynamo. 

In what follows, we restrict ourselves to the bipolar 
topology and assume that this field, although it has a 
profound dynamical influence on the disc, is not strong 
enough to significantly perturb its radial balance. Thus, 
the disc is supposed to be geometrically thin (its half- 
width at a distance r verifies h{r) <C r), in a quasi Kep- 
lerian equilibrium. The magnetohydrodynamic equations 
describing stationary, axisymmetric MAES are then the 
following: 



V • pit = 

pu ■ Vu ~ —pV^G 



VP + J X B + V -T 



(1) 

(2) 
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J = —V X B 



(3) 
(4) 
(5) 



where $g = —GM/{r^ + z^)^/^ is the gravitational po- 
tential of the central object (disc self-gravity is neglected) 
and T is a viscous stress tensor (Shakura & Sunyaev 
1973). Both disc velocity and magnetic field were de- 



composed into poloidal and toroidal components, namely 
u = Up + flre^ and B — Bp + B^e^ respectively. A mag- 
netic bipolar topology imposes an odd B^ with respect to 
z, as well as a poloidal field 



Bp — -Va X Ba 

r 



(6) 



where a(r, z) is an even function of z. A magnetic surface, 
which is a surface of constant magnetic flux, is directly la- 
belled by a(r, z) = ao{ro, 0). The distribution of magnetic 
flux through the disc is a free (and unknown) function. 
We will therefore use a prescription consistent with the 
complete set of equations. 

Since the magnetic field threads the disc, steady-state 
accretion requires that matter diffuses through the field. 
As usual in astrophysics, normal transport coefficients are 
far too small to account for the expected motions. Drift 
between ions and neutrals, known as ambipolar diffusion, 
could play such a role in accretion discs (Konigl 1989| , 
Wardle & Konigl 1993] ). However, jets are observed in 
a wide variety of objects, which suggests that they are 
produced by a mechanism independent of the disc ionisa- 
tion degree. Moreover, it is now well known that magne- 
tized discs are prone to instabilities (e.g. Balbus & Hawley 
1991| , Tag ger et al. |1992| , Fo glizzo & Tagger |1994 Curry 



& Pudritz 1995, Spruit et al. 1995). Thus, we assume that 
the disc is turbulent and that the non-linear evolution of 
this turbulence provides the required anomalous transport 
coefficients, such as the magnetic diffusivity i>m, resistiv- 
ity rim = iJ'oVm and viscosity Vv (appearing in T). All 
these transport coef ficien ts should achieve a comparable 
level (Pouquet et al. 1976|) . Nevertheless, we introduced in 
Eq.(5) a "toroidal" magnetic diffusivity v'^ to account for 
anisotropy with respect to the "poloidal" diffusivity 
Indeed, the winding up of the field due to disc differential 
rotation can lead to strong instabilities and thus, to en- 
hanced transport coefficients with respect to the toroidal 
field (see FP95). 

Finally, in order to close the system, we use a poly- 
tropic approximation 



P = JCp' 



(7) 



with /C constant along a magnetic surface, in the isother- 
mal case (7 = 1). Magnetized jets driven by an accretion 
disc can be initially launched either by predominant mag- 
netic effects (magnetically-driven) or by enthalpy alone 



(thermally-driven). For discs without a hot corona (cold 
MAES with negfigible enthalpy), FP95 showed that jets 
could indeed be produced through magnetic effects alone. 
For these magnetically-driven jets, the vertical profil of 
the temperature has not a strong influence on plasma dy- 
namics, therefore allowing a description with an isother- 
mal structure. The reason lies in the fact that most of the 
accretion power goes into the jets as an MHD Poynting 
flux. Thus, magnetically-driven jets will always be associ- 
ated with weakly dissipative discs, where thermal effects 
can be neglected (see below). 

Such a cold MAES corresponds to a "clean" magnetic 
structure at the disc surface, with just one polarity. In fact, 
it is likely that the magnetic structure would consist of 
open field lines coexisting with small scale loops anchored 
at different radii (Galeev et al. 197£, Heyvaerts & Priest 



1989). Such a situation would then most probably result 



in the formation of a hot corona. The treatment of such a 
hot MAES is postponed to future work. 

The local state of a cold MAES is mainly characterized 
by the set of following parameters evaluated at the disc 
midplane, 



r 

rur 



VAh 
B^ 



A = 



d In Ma 
dlnr 
(J X B)^ 





(8) 



(V-T) 



2gM 



namely, the disc aspect ratio e ^ 1, the magnetic 
Reynolds number TZm, the turbulence level parameter a™, 
the disc magnetization ^, the ejection index ^ > and the 
ratio A of the magnetic to the viscous torque (FP95). This 
ratio depends on the magnitude fj, of the field, the mag- 
netic shear q of order unity, defined as 



h dB. 



Bo dz 



(9) 



2=0 



and being a measure of the toroidal field at the disc sur- 
face (Ferreira & Pelletier 19934 hereafter FP93a) and, of 
course, on the magnitude of the viscosity through the well 
known parameter (Shakura & Sunyaev 1973). Under 
the sole assumption that the disc is in Keplerian balance, 
angular momentum conservation implies that 



1 -I- A ~ 7^„ 



(10) 



must be verified at the disc midplane. If large scale mag- 
netic fields are irrelevant (A <C 1), one recovers that the 
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Reynolds number TZy — rur/vv is of order unity, which is 
required in standard viscous discs. On the contrary, when 
they significantely brake the disc, one obtains that the de- 
gree of curvature of the field lines at the disc surface (mea- 
sured by TZm) depends mostly on how strong is A. BP82 
showed that in order to magnetically launch jets without 
a hot corona, the magnetic structure needs to be bent 
by more than 30° with respect to the vertical axis. This 
implies a magnetic Reynolds number of order TZm ^ £^^, 
therefore A <; r/h. Cold jets require then a dominant mag- 
netic torque, imposing a corresponding value on both the 
field strength and shear. 

Energy conservation equation (FP93a), 



Pi 



lib 



2P 



2P, 



MHD 



2Pf 



th 



(11) 



shows that the total available mechanical power Pm, is 
shared by the disc luminosity Prad, the outward MHD 
Poynting flux Pmhd and the flux of thermal energy Pth 
from each surface of the disc. The liberated power is suc h 
that Pub <, Pace = GMMae/^n (scc FP95, Ferreira |l996|) , 
where is the disc inner radius and Mae is the accretion 
rate at the disc outer edge Vg. Around a compact object, 
this accretion power is 



p 

ae 



4.7 10 



45 




(12) 



— 3rg, where = 2GM is the 



for Af, = lMQyr~\ r 
Schwarzschild radius. Around a protostar, this power is 



Pa. 



10 



32 



Mae 



M 



.5M, 



— ) erg s 



(13) 



for A/, = lO~^Af0yr"^, — 30rQ (10 stellar radii for a 
typical T-Tauri star). Since the available energy is stored 
into rotation, energy conservation can be rewritten as 



2-frad 



2R 



jet 



1-e 
TTa 
A+_e 
TTX 



Pub 
Pm 



(14) 



where O < 1 is the fraction of the internal energy that 
was transferred in the jet as enthalpy. The total jet power, 
which can be either radiated away or injected at the ter- 
minal shock (e.g. extragalactic radio lobes), arises from 
the sum of the MHD Poynting flux and the plasma ther- 
mal power. Thus, magnetized accretion discs without a 
hot corona (A ~ r/h, = 0) radiate only a fraction of 
order h/r of the jet power. 

2.2. From accretion to ejection 

MAES are intricate structures where accretion and ejec- 
tion are interdependent. Therefore, we will frequently 



make references to disc physics and quantities in our in- 
vestigation of jet physics. Let us then, for the sake of com- 
pleteness, summarize here results concerning magnetized 
disc physics. 

Two simultaneous processes are responsible for accre- 
tion (FP93a): (1) a turbulent magnetic diffusivity al- 
lowing matter to steadily diffuse through the field; (2) 
a dominant magnetic torque Ffj, = J^Br — JrBz, which 
brakes the disc and stores into the field both angular mo- 
mentum and mechanical energy of the plasma. A local 
quasi-magnetohydrostatic vertical equilibrium is achieved, 
the plasma pressure gradient hardly competing with both 
tidal compression and vertical magnetic pinching force. A 
radial magnetic tension slightly counteracts gravity, thus 
leading to a sub-Keplerian rotation rate. 

As a result, plasma inside the disc (see Fig. |]) is be- 
ing accreted, slightly converging towards the disc mid- 
plane (both Ur and Uz are negative). How then is ejection 
achieved ? 

Ejection comes out naturally if the radial current den- 
sity Jr decreases vertically on a disc scale height (Ferreira 
& Pelletier [l993b| ). Indeed, in a Keplerian disc, this is the 
only possibility to change the sign of the magnetic torque 
F^. This torque must become positive at the disc surface 
in order to provide a magnetic acceleration. Such a sit- 
uation requires that the counter current due to the disc 
differential rotation balances the current induced by the 
unipolar induction effect (FP95). If this is not fulfilled, the 
complete structure is unsteady: in fact, the requirement of 
stationarity provides the level of the "toroidal" diffusivity 

The vertical decrease of Jr leads to the following im- 
portant effects: 

(1) the vertical Lorcntz force decreases, allowing the 
plasma pressure gradient to gently lift up matter from the 
disc surface; 

(2) matter is azimuthally accelerated by the magnetic 
field, leading then to a radial centrifugal acceleration. 
Note that, although it is plasma pressure that provides a 
positive vertical velocity, this process is purely magnetic. 
It arises naturally if the magnetic pinching force (that 
must be comparable to the plasma pressure gradient) de- 
creases vertically. 

In order to get an insight on the jet physical processes, 
it is worthwhile to project the Lorcntz force parallel and 
perpendicular to any poloidal magnetic surface. 





JDri 








27rr " 








27rr 







(15) 



27rr 



-V,/ 



Here, / = 2'KrB^/^o < is the total current flowing 
within this magnetic surface, V|| = {Bp ■ V)/Bp and 
= (Va • V)/|Va|. This shows directly that plasma 



Jonathan Ferreira: Magnetically-driven jets from Keplerian accretion discs 5 




Fig. 1. Side view of a magnetized accretion disc driving jets, with e — 10" , a,„ = 1 and an ejection efficiency ^ — .004. Plasma 
(dotted lines) enters the structure at its outer edge and is accreted with a slight converging motion through the magnetic field 
(solid lines). Since the magnetic pinching force decreases vertically, plasma reaches a layer where the pressure gradient slowly 
expells it. Diffusion is still necessary at the disc surface in order to allow a transition between the accretion disc and the ideal 
MHD jet. This transition region is approximatively one disc scale height thick. Farther out, plasma is frozen in a particular 
magnetic field line and is accelerated through the Lorentz force. Note the scaling factor applied to the vertical axis. 



is accelerated by the current leakage through this surface 
(V||/ > 0). This effect gives rise to both a poloidal (F||) 
and a toroidal (F^) magnetic force, the latter providing 
the centrifugal force. Thus, cold jets are better referred 
to as being magnetically-driven rather than centrifugally- 
driven (this was also pointed out by Contopoulos & 
Lovelace 1994 ). When the current / vanishes, or when it 
flows parallel to the magnetic surface, no magnetic accel- 
eration arises anymore and the plasma reaches an asymp- 
totic state. Note also that the way this current is dis- 
tributed across the magnetic surfaces is of great impor- 
tance for the jet transverse equilibrium {F±). This shows 
that one has to be careful when dealing with a particular 
current distribution, since it is of so great importance in 
both jet collimation and acceleration (Appl & Camenzind 



1993a, 1993b). 



The global current topology is linked to the ejection ef- 
ficiency, which is measured by the ejection index ^ = 
in a viscous disc without jet). For small ejection efficien- 
cies (^ < 1/2), ejection takes place against a pinching 
vertical Lorentz force, with a positive radial component. 
For those tenuous ejections, the current enters the disc 
at its inner edge, flows up inside the jet and closes back 



along the axis. High ejection efficiencies (^ > 1/2) are 
achieved when matter is lifted by both plasma and mag- 
netic pressure gradients (positive vertical Lorentz force), 
the corresponding radial Lorentz force being negative. In 
this case, the current flows down the jet, enters the disc at 
its surface and closes in an outer cocoon (Ferreira 1994). 

The set of MAES parameters (8) is diminished by 
the requirement that the overall structure is in steady- 
state. This implies that the flow has to smoothly cross 
the usual MHD critical points it encounters (see FP95). 
The constraint due to the first critical point, the slow- 
magnetosonic (SM) point, allows these two current topolo- 
gies to be realized. In what follows, we will show that only 
tenuous ejections can produce trans- Alfvenic jets. 



3. Non-relativistic, magnetically-driven jets from 
Keplerian discs 

3.1. Governing equations 

Above the disc, both the decay of the turbulent diffusivity 
and the action of the poloidal Lore ntz f orce leads mat- 
ter into ideal MHD state (Ferreira 1996 ). In this state, 
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matter is frozen in the field, the poloidal motion occur- 
ring along a magnetic surface. Together with axisymme- 
try, this allows to describe jets as a bunch of magnetic 
surfaces, nested on the accretion disc at an anchoring ra- 
dius Tq. In this section, we write the well-known govern- 
ing equations for self-confin ed jet s in a form similar to 
the one Pelletier & Pudritz ( |1992| , hereafter PP92) used. 
Written in such a form, these equations will allow us to de- 
rive general properties and conditions for non-relativistic, 
magnetically-driven jets when one takes into account the 
underlying disc. In Sect. 4, we will come back to the set 
of equations as described in Sect. 2 and solve them from 
the disc equatorial plane to the "jet end" . 
In the ideal MHD jet, Eq.(4) becomes 



Ur, 



via) 



B 



p ' 



(16) 



where ri{a) is a constant along a particular magnetic sur- 
face. Its value, T] = y/HoPA, is directly evaluated at the 
Alfven point (labelled with the cylindrical coordinates ta 
and za), where the jet poloidal velocity reaches the local 
Alfven speed. Equation (5) becomes 



r2*(a) = — r/' 



' P'opr 



(17) 



and defines the rotation rate 11* of a magnetic surface. 
Since the field lines are anchored on the accretion disc, 
they rotate with roughly the same rotation rate as matter, 
namely 51* ~ f2o. Thus, one can look at this equation as 
providing the amount of toroidal field required to maintain 
matter, which is frozen in the field, rotating in the jet at 
a different rate than the field. 

Using these two equations, the steady-state angular 
momentum conservation equation can then be written as 



rB^ 



(18) 



where ri*r^ is the total specific angular momentum car- 
ried away by both matter and field in a particular mag- 
netic surface. The radial and vertical momentum con- 
servation equations are usually replaced by their projec- 
tion along (Bernoulli equation) and perpendicular (Grad- 
Shafranov equation, also known as transfield equation, 
Tsinganos |1981 ) to such a magnetic surface. For poly- 
tropic flows, Bernoulli equation reads 



u 

y 



H 



- a 



rB^ 



E{a) 



(19) 



where E{a) is the constant specific energy carried by the 
jet and the jet enthalpy H is defined as \/H = VP/ p. For 
magnetically-driven jets (i.e. cold), this enthalpy plays no 
role neither in jet formation nor in its coUimation. Indeed, 
it is only a fraction of order of the Keplerian speed 
and so it will be simply dropped in the following section 



(but note that the solutions in Sect. 4 are obtained with 
the full set of self-similar MHD equations). As shown in 
FP95, ejection of matter from the disc arises naturally 
due to the vertical decrease, on a disc scale height, of the 
radial current density. This is a pure magnetic process, 
without any help from thermal effects. On the contrary, 
thermally-driven jets require a high enthalpy with a poly- 
tropic index 7 variable along a magnetic surface, namely 
close to unity near the disc to account for coronal heat- 
ing and clos er to the adiabatic value furt her a way (We- 
ber & Davis 1967 , Tsinganos & Trussoni 1991 , Sauty & 
Tsinganos 1994). This last possibility is postponed to fu- 
ture work. Bernoulli equation (19) describes how the total 
energy carried by the flow is transformed into kinetic en- 
ergy, for a given magnetic conflguration. Hence, one can 
interprete this equation as providing the velocity matter 
reaches in a given "magnetic funnel" . The shape of this 
funnel, or more precisely the jet transverse equilibrium, is 
provided by Grad-Shafranov equation 



V • (rn^ - 1) = p _ - r!— ^ 

/ior^ 1^ da da 



dil* 
da 

Bl + m'BldlnpA 
2/io da 



(20) 



where we introduced the Alfvcnic Mach number = 

The set of equations (n6|) to (20), together with Eq.(R), 
completely describe magnetic jets where plasma pressure 
plays no dynamical role. Hereafter, we describe the leading 
parameters of these equations. 



3.2. An unique parameter for cold jets 

We choose to characterize our jet solutions by using pa- 
rameters as defined in the pioneering work of BP82. Since 
each parameter is strictly defined at the footpoint Tq of a 
magnetic surface, it characterizes the jet state only locally. 
Thus, in a realistic 2-D case, one would have to prescribe 
them for a range of anchoring radii in the disc. However, 
if we assume that jets are produced from a large radial ex- 
tension in the disc, then we can look for a jet solution de- 
fined with parameters that are constants (as in self-similar 
solutions), or slowly varying with the radius. Through all 
this paper, we will label with a subscript "o" any quan- 
tity evaluated at the disc midplane and a subscript "A" 
at the Alfven surface. In particular, ~ p{ro, 0) will be 
the midplane density and Uo = — Ur(ro,0) the accretion 
velocity at a radius rg. 

The first parameter describes the magnetic configura- 
tion, namely 



/3 = 



din a 
din To 



(21) 
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with < (3 < 2 (FP93a). Jets with constant /3 have a mag- 
netic field varying with a power law of the radius. BP82's 
self-similar solutions were obtained with the prescription 
(3 — 3/4:. The second parameter, defined as 



A 



(22) 



is a measure of the magnetic lever arm that brakes the un- 
derlying accretion disc. Since the available energy is stored 
in the disc as rotational energy, a constraint on this lever 
arm arises through the requirement that jets are powered 
with a positive energy. Indeed, the total specific energy 
carried away by the (cold) magnetic structure writes 



E{a) = £{a) + nlr\ 



(23) 



where £{a) ~ 
tial well if the 



jT^/'i,. Jets become free from the poten- 
"rotator" energy, namely Q,^r\, overcomes 
the generalized pressure £ (PP92). This is fulfilled pro- 
vided A > 3/2. The last parameter. 



s7 



(24) 



measures the mass load on a particular magnetic surface. 
Since cold jets carry away the whole disc angular momen- 
tum, one expects to find a systematic relation between 
mass load and lever arm. 

These parameters must be linked to the ejection index, 
since it is a local measure of the ejection efficiency. A way 
to find such a link is to look at the ratio of the MHD 
Poynting fiux to the kinetic energy flux along a magnetic 
surface. 



(25) 



Such a ratio measures the amount of energy stored as 
magnetic energy with respect to the kinetic energy, being 
therefore a measure along the jet of the acceleration effi- 
ciency. At the jet basis, identified here as the SM-surface, 
the disc provides 



crsM 



Ma 



(26) 



where we used mass {dMa/dr — 2dMj/dr — Airpuzr) and 
angular momentum conservation. Thus, the ejection index 
is also quite accurately the injection of energy into the jets. 
While the magnetic lever arm is directly A = 1 -I- (T5j\,//2, 
the mass load writes 



K — 2(7 



SM 



B. 



4>,SM 



Bo 



(27) 



which provides k ~ ^ ^ (see Appendix Thus, 
the leading parameters for magnetically-driven jets can 



be simply expressed as 



A 



1 



(28) 



The expression for the magnetic configuration parameter 
(3 is general for any Keplerian disc, allowing to take into 
account simultaneously magnetic and gravitational terms 
(FP93a, FP95). All three parameters for cold jets can be 
quite accurately replaced by ^: magnetic configuration and 
both mass load and lever arm are therefore tightly related. 
As a consequence, the parameter space for cold jets is 
mostly controlled by the range of allowed ejection indices. 

By using (28), we are now able to see that any jet 
model that would allow ejection with a high efficiency 
(namely, ^ > 1) would not obtain free jets {E{a) < 0). 
Nevertheless, jet solutions that successfully cross the SM- 
point with ^ > 1 were found to be possible (Ferreira 1994). 
Such a situation is therefore a transient feature, matter 
failing to reach the Alfven surface and falling down to the 
disc after having been ejected out. 

3.3. Constraints on the ejection index 

3.3.1. Minimum ejection index: disc vertical equilibrium 

A look at the parameter space in FP95, derived with the 
sole constraint due to the slow-magnetosonic point, indi- 
cates that there is a minimum ejection index ^min- Here, 
we show that this minimum ejection efficiency is related to 
the existence of a quasi-magnetohydrostatic (MHS) equi- 
librium inside the disc. Such an equilibrium, described by 



du^ aP _ ^ _ 9 Bf.+ B^ 

^ dz dz ^ dz 2po 



(29) 



is obtained when plasma pressure gradient balances both 
tidal compression and magnetic squeezing. This magnetic 
squeezing depends on the value of the magnetic field, mea- 
sured by the parameter fj, (see Eq.(8)), and on both curva- 
ture (for Br) and shear (for B^) effects. When the ejection 
efficiency decreases, the lever arm increases and so does 
the magnetic compression due to field curvature. The only 
way to find out another equilibrium state is then to de- 
crease the value of the field itself (that is, n). However, 
this results in an enhancement of the magnetic squeezing 
due to shear. 

This can be qualitatively shown like follows (for more 
details, see Appendices ^ and |^). The induction equation 
(5) can be written in the disc as 



v'mJr — VoJo + r dzBp ■ Vf2 - B^Uz 



(30) 
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where the last term describes the effect of advection. In the 
resistive disc, this term is negligible in comparison with 
the others. Thus, the toroidal field at the disc sm'face. 



B 



(31) 



is mainly the sum of two contributions, B^" — —fj,oJoh 
due to the unipolar induction effect and Bj^ , the counter 
current (of positive sign) provided by the disc differential 
rotation. If the disc were rigidly rotating (i.e., a Barlow 
Wheel), it would generate the first contribution but not 
the second one, therefore producing no jet. When fj, de- 
creases, the magnitude of this counter current decreases, 
thereby increasing the toroidal field at the disc surface and 
with it, the pinching effect of the corresponding magnetic 
pressure gradient. 

There is thus a S,min below which no quasi-MHS equi- 
librium can be found anymore. The exact value of this 
minimum ejection index is impossible to find analytically, 
since it depends on a subtle equilibrium between terms 
that are all of the same order of magnitude. When such 
equilibrium is crudely treated, one can in principle still 
obtain a matching with jet solutions but the disc would 
certainly not survive an overwhelming magnetic squeez- 
ing. 

In Section 4, we show the parameter space obtained for 
self-similar solutions, thus providing a numerical value for 
Cmm- Because self-similarity does not influence solutions 
close to the Keplerian disc but allows instead to take into 
account all dynamical terms, we believe that this value 
is general. Since what limits the minimum ejection index 
is the increasing toroidal field, one can only significantly 
lower ^min by decreasing the value of B^° , that is, by 
lowering Jq. Therefore, only discs where both viscosity 
and magnetic torques are relevant (A ^ 1, requiring a 
hot corona to produce jets), should allow ejection indices 
much smaller than those displayed here. 

3.3.2. Maximum ejection index: jet acceleration 

Conservation of angular momentum (Eq.([l^)) implies that 
once a large scale magnetic field steadily brakes the disc 
with a significant torque, matter must be accelerated up 
to the Alfven point. This severily constrains the mass load 
in the jet, thus providing an upper limit on the ejection in- 
dex. When gravity becomes negligible, Bernoulli equation 
can be rewritten (FP93a) as 



It measures the discrepancy between the angular veloc- 
ities of matter and magnetic surface on which matter 
flows (PP92). At the jet basis, its value is approximately 
zero, then it increases and reaches unity for highly super- 
Alfvenic jets. The Alfvenic Mach number becomes then 
simply 



,Bo 

Br, 



(34) 



A necessary condition for trans-Alfvenic jets is obtained 
by requiring that the above expression reaches or becomes 
much bigger than unity. At this stage, it is convenient to 
introduce the ratio, measured at the Alfven point, of the 
rotation velocity of t he ma gnetic surface to the poloidal 
Alfven speed (Michel |l969| , PP92) 



LUA 



Vap,a 



Br, 



B. 



p.A 



(35) 



This fastness parameter is a useful quantity since it en- 
closes an information that we cannot have without solving 
first the whole structure, namely the angle of the poloidal 
magnetic field with respect to the vertical axis. This ratio, 
which must verify 



UJA > 



B, 



Br, 



(36) 



will allow us to derive a general condition for trans- 
Alfvenic jets, in two extreme cases. 

For "powerful" jets (high (Jsm)^ matter is expelled off 
the disc carrying almost no angular momentum. Accelera- 
tion takes place only if a a < o'sm, with a a 
these jets require 



2uj\. Thus, 



1 

2? 



(37) 



This corresponds to the extreme case where it takes almost 
no power from the magnetic structure to accelerate plasma 
up to the Alfven surface. Such a configuration could be 
obtained in the low mass load limit. 

For "weak" jets (low asAi), the Alfvenic Mach number 
has to reach at least unity. At the Alfven point, Bernoulli 
equation writes a\ — (1— .g^)A— 3. This shows that matter 
reaches its maximum velocity when gA = 0, namely when 
the toroidal field is almost zero. Using Eq.(^), it can be 
seen that these jets require 



(32) 



where the function g, defined as O = ri*(l — g), can be 
expressed as 



(33) 



— 1 



X 



X 



4 A - 3 



(38) 



This general condition shows that magnetically-driven jets 
are fast rotators {wa > 1) in order to successfully reac h 
the Alfven surface (BP82, PP92, Rosso & Pelletier |1994D . 
This is required even in this case, where matter barely 
reaches it despite the almost complete transfer of energy 
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from the magnetic structure. High mass loads could give 
rise to such jets. 

Compiling criteria ( p7| ) and one gets the neces- 

sary condition for trans- Alfvenic jets 



1 + 2C 2 1 



(39) 



which must be satisfied at each magnetic surface. Here, 
we get a much more precise constraint since it requires 



(40) 



This implies that higher ejection efficiencies are com- 
pletely inconsistent with a steady trans- Alfvenic regime. 
Note that this value should be viewed as an upper limit 
for ^, the maximum ejection index S^max being in fact 
smaller. This is due to the weak constraint we used in 
the high-asM case. Contopoulos ( |1995D proposed a sim- 
ple jet model where the magnetic field has only a toroidal 
component. In such a case (obtained in the limit k very 
large), the Alfven singularity disappears and the driving 
mechanism for jet acceleration is the magnetic pressure 
gradient. We would like to stress that such a configura- 
tion is, in principle, possible as a transient outburst from 
thin discs. Indeed, FP95 showed that magnetically-driven 
jets can be viewed as being either "centrifugally-driven" 
for small ejection indices (^ < 1/2, the vertical compo- 
nent of the Lorentz force pinching the disc) or "magnetic 
pressure-driven" for high ejection indices (^ > 1/2, the 
toroidal magnetic pressure lifting matter up, consistent 
with the large k limit used). However, only the first sit- 
uation allows a steady state with respect to the Alfven 
surface. The possibility remains that such a configuration 
would be steadily produced by a thick disc. In the next sec- 
tion, we compute global self-similar solutions, derive both 
£,min and (,max numerically and obtain the full parameter 
space for cold MAES. 

4. Trans- Alfvenic, global disc-jet solutions 

The self-similar ansatz 

In order to construct global solutions from the disc equa- 
torial plane up to a jet asymptotic regime, the full set of 
non-linear, partial differential equations (1) to (5) along 
with Eq.(|) must be solved. Since MAES are intrinsically 
2-D, this requires either a numerical approach (not yet 
available) or some method to reduce this set to a set of or- 
dinary differential equations (but enforcing some symme- 
try to solutions). This is what allows self-similarity, which 
is the reason why it has been used by many authors (see 
Tsinganos et al. ( 1996 ) for a review). 

We look then for self-similar solutions so that any 
quantity Q can be written 



Q(r,z) 



/q(^) 



where is the external radius of the disc (a standard vis- 
cous disc is probably established from r,, to larger radii). 
The self-similar variable is chosen to be a; = z/h{r), con- 
sistent with the gravity field at the disc neighbourhood, 
but also everywhere since h{r) = er. Such a form, h oc r, 
is expected to arise in a Keplerian disc pervaded by a 
large scale magnetic field (FP93a, FP95). The magnetic 
flux function is chosen to be 



a(r, z) 



13 



(42) 



where ij}{—x) = ip{x) and the pressure prescription such 
that JC[a) =/Ce(ro/re)"'^, with a^; = (1 - 7)(2/3 - 3) - 1, 
the polytropic index 7 remaining a free parameter. Since 
we focus on cold jets, its value has no importance and we 
choose 7 = 1 for convenience (see discussion after Eq.(0)). 

The set of PDEs can be separated into two sets, the 
first consisting of algebraic equations between the indices 
aq and /3, the second of non-linear ODEs on the functions 
/q, V and their derivatives (see Appendix B in FP95). 
Within such a prescription, all the dynamical terms can 
be included, allowing henceforth a complete study of the 
physics of ejection from Keplerian discs. However, only 
non-relativistic jets can be studied within the same self- 
similarity prescription. Indeed, relativistic speeds enforce 
a scaling with the speed of light (Li et al. 1992| ), which is 
inconsistent with disc physics. 

Since the boundary values /q(0) at the disc midplane 
are known, we just have to integrate the set of ODEs from 
a; = to infinity. Once we have a global solution as a 
function of a: = z/er, it is straightforward to obtain the 
variation of any quantity along any magnetic surface a = 
constant with 



Q{x{a)) ^ QofixWx)-"'^/^ 



(43) 



where Qo — Qe(?'o/''e)"'' (note the use of the notation 
x(a) in this case). The main difficulties come from the 
presence of three critical points, whenever the plasma ve- 
locity V = u ■ n, where 



(44) 



equals the three usual MHD phase speeds: slow- 
magnetosonic, Alfven and fast-magnetosonic, in the di- 
rection n. Self-similarity not only modifies the definitions 
of these phase speeds, but constrains also the direction 
of propagation where critical points appear (see FP95 for 
more precisions). Hence, the critical velocity V is roughly 
Uz at the disc neighbourhood where the flow becomes 
super-SM, and Ur much farther away. 



^^2^ -J 4-2. MAES parameter space 
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Fig. 2. MAES parameter space for an isothermal structure 
with Om = 1, for various disc aspect ratios: e = 10~^ (solid 
line), 10"^ (dotted line), 10"^ (short-dashed line) and 7 10"* 
(long-dashed line). The disc magnetization jj, and the magnetic 
Reynolds number TZrn were obtained as regularity conditions, 
in order to get respectively trans-SM and trans-Alfvenic so- 
lutions. The other two parameters, the magnetic shear at the 
disc surface q and fastness parameter ua are calculated. As 
^ decreases, fi must decrease in order to balance the increase 
in magnetic shear q and curvature TZm£ (see Sect. 3.3.1). It is 
noteworthy that, although the disc parameters change a little, 
the main effect of decreasing e is to shift the range of allowed 
^ to higher values. We confirm here that trans-Alfvenic jets 



require ua^I (PP92, Rosso & Pelletier 1994) 



4.2.1. The slow-magnetosonic point 

We obtain trans-SM solutions by adjusting the strength 
of the field for a given set of parameters (FP95). Thus, for 
fi bigger than a critical value fic, the magnetic compres- 
sion is too strong and leads to a vanishing density (hence 
infinite vertical acceleration). When /i < /Zc, the magnetic 
compression is insufficient and too much mass is expelled 
off, which eventually falls down again. 

Thus, the smooth crossing of the first critical point is 
directly related to the disc vertical balance. Following Li 



(1995), we can express that a necessary (but not sufficient) 
condition for stationarity is 



where the subscript refers to the disc surface. In terms 
of MAES parameters, this condition for cold jets becomes 



< 



47^i£2 



where the function fj^ is of order unity (see Appendix 
This implies fi < 3/2, which shows that jets can- 
not be produced from magnetically dominated discs. Since 
plasma pressure plays such an important role in sustaining 
the disc, only structures close to equipartition are steady 
(FP95). 

The minimum ejection index is found by re- 

quiring that all solutions become super-SM. We found 
Cmin = 0.004 for an isothermal disc with e = 0.1 and 
a^n = 1. If the disc aspect ratio e (and also the turbu- 
lence level am, see below) decreases, the influence of ad- 
vection in Eq.(|o|) also decreases and the radial current 
density Jr is smaller, hence providing a smaller value of 
the toroidal field at the disc surface (/+, see Appendix 
As an example, we obtain /+ ~ 0.6,0.5,0.45 for 
e = 0.1,0.01,0.001 respectively. The magnetic compres- 
sion decreasing, the range of allowed ejection efficiencies 
is shifted to higher values (see Fig. ||). Moreover, by using 
an isothermal prescription for the disc temperature, we 
have underestimated the "lifting" efficiency of the plasma 
pressure gradient. Thus, ^min = 0.004 should be seen as 
the lower limit for cold jets from Keplerian discs. However, 
since ^ < 0.15 is a general constraint, one can expect to 
find that below a certain value of e (and a™), no steady- 
state solutions can be found anymore. 

The work presented here strongly differs from the work 
done by Wardle & Konigl ( |1993D and Li ( |1995D . These au- 
thors found no limiting ^ and were thus able to obtain 
jet solutions with an enormous lever arm A, and an ar- 
bitrarily small mass load k. Here, it has been found that 
such a situation is impossible to achieve in steady-state. 
This disagreement comes from their different treatment 
of the highly sensitive disc vertical equilibrium. Indeed, 
Wardle & Konigl investigated the vertical structure at a 
given radius and made the assumption (valid for the jet) 
pUz = constant in order to deal with the mass conser- 
vation equation. Li used a self-similar approach, but did 
not make a smooth transition between the disc and its 
jets, crudely matching the disc MHS solutions to jet ones. 
He however found that the mass load was very sensitive 
to the plasma rotation rate at the disc midplane Qq: the 
smaller flo the smaller k. This can now be understood with 
what was previously said. Indeed, flo = ilfc(l — lUq) with 
uJo — fJ.TZms'^ f^f = o(e). Thus, a smaller fio is achieved 
with a bigger e, which implies a smaller ejection efficiency 
(hence, k ~ ^). 



(45) 4.2.2. The Alfvcn point 



(46) 



In order to get trans-Alfvenic solutions, we adjust the 
magnetic Reynolds number TZm-, which controls the bend- 
ing of the poloidal field lines at the disc surface. For TZm 
bigger than a critical value TZc, the bending is too strong 
and leads to an overwhelming centrifugal effect, such that 
rir^ > r2,r^ (i.e., becomes unphysically positive). On 
the contrary, if TZm < Ti-c, the centrifugal effect cannot bal- 
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ance the magnetic tension, which leads to an unphysical 
closure of the magnetic surfaces {Br becomes negative). At 
each trial for TZm one has to find another critical that 
allows a super-SM solution. Note that fixing TZm is the 
same as fixing the accretion velocity at the disc midplane, 
for any given magnetic diffusivity. We have the freedom to 
do it because this accretion velocity is also controlled by 
the magnetic field. If we were in a situation where accre- 
tion is mainly due to the viscous stress, then TZ^n would 
also be given. 




Fig. 3. Parameter space for cold jets launched from an isother- 
mal accretion disc, with am = 1 and for various disc as- 
pect ratios: e = 10~^ (solid line), 10~^ (dotted line), 10"'^ 
(short-dashed line) and 7 10~^ (long-dashed line). Due to the 
SM constraint, we get here a smaller parameter space as in 
BP82. As showed, the jet parameters do not significantly vary 
with any parameter but ^. 



Figure ^ shows the parameter space for cold jets in 
the K-X plane for am = 1- Solutions with k > 0.1 do 
not allow super-Alfvenic jets (BP82, Wardle & Konigl 
1993). We found here that ^max = -08, the maximum 



Alfvenic Mach number reached by the flow being barely 
bigger than unity. For disc aspect ratios below 5 10""*, 
no steady-state solutions are found, £,min becoming equal 
to (,max- Thus, trans-Alfvenic jets are produced for discs 
where 5 lO'^ < £ < 10-\ the upper bound arising from 
the thin-disc approximation. 

Using Eq.(p^, we can express the density at the Alfvcn 
point as a function of the disc midplane density and MAES 
parameters, namely 



Po 4 



(47) 



This relation is local and independent of any boundary 
condition that could constrain the jet behaviour. Since 



we approximately obtain numerically such a value for the 
Alfven density, as well as a consistency with all our general 
requirements of Sect. 3.3.2, we believe that self-similarity 
doesn't affect the solutions up to this point. Therefore, 
our parameter space should be viewed as general for cold 
jets launched from Keplerian discs. 




o □ 
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Fig. 4. Low-efficiency (^ = 0.005, left pannels) and high-effi- 
ciency (f = 0.05, right pannels) trans-Alfvenic jets for 
e = 10~^, am = 1. The upper pannels show the angular mo- 
mentum transfer from the field (>, dashed line), to the plasma 
{jm, solid line). At the disc surface (x = 1), all the specific 
angular momentum Q,r\ is carried by the field, but it is af- 
terwards completely transferred into the plasma. Below, the 
corresponding function g = 1 — fJ/il, is displayed for both 
cases. The cross labelled with "A" tells the position of the 
Alfven surface. In the low-efficiency limit, the jet reaches it 
with p ~ 1 {uA ~ 1.8), most of the power being still carried by 
the field. In the high-efficiency case, the flow requires already 
almost all the power to reach Alven speeds {coa = 1.19). 



At this stage, one cannot tell with precision what will 
be the asymptotic state of the jet (unless of course by 
propagating the solution, as we will do next section) . How- 
ever, by using general arguments based on the work of 
Heyvaerts & Norman ( 1989| ), one can give clues on the 
degree of collimation achieved, according to the amount 
of current I a 



Ia 
Ism 



9a{1 



9A 



(48) 



still available at the Alfven point. Hereafter, we make 
a distinction between two kinds of jets, depending upon 
their state at the Alfven surface (see Fig. ^). 

Current-carrying jets {gA ^1) reach the Alfven sur- 
face while most of the angular momentum is still stored 
in the magnetic structure. They are powerful (correspond- 
ing here to "low" ejection indices 0.004 < ^ ^ 0.025) and 
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could reach highly super- Alfvcnic speeds (m^ ^1), with 
n ~ ri*r^/r^ and a maximum poloidal velocity 



(49) 



Such jets could allow a cylindrical coUimation if this cur- 
rent does not vanish completely. 

Current-free jets {gA ^ 1) become super-Alfvenic at 
the expense of almost all the available current. They cor- 
respond to "high" ejection efficiencies 0.025 ^ ^ < 0.08 
(for am = 1, see below). According to Heyvaerts & 
Norman, only parabolic coUimation could be asymptot- 
ically achieved. These jets could reach moderate speeds 
(m^ > 1), with ^ and a maximum poloidal velocity 



r!or„(A- 3)1/2 



(50) 
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Equatio n (|35 ) provides, for small ejection indices (Pel- 
letier et al. |1996), 



-1/2^-1/2 



(51) 



The fastness parameter lua strongly depends on both the 
shear parameter q (henceforth on am) and the ejection 
index ^. The highest value of uja (hence, of qa) will be 
achieved with am — 1 and the smallest value of ^. Since 
the disc vertical equilibrium limits the latter, lua is there- 
fore also hmited {1 < loa < 2, see Figs. ^ and ||). 

As am decreases, the influence of plasma on the field 
grows: the counter current due to the disc differential ro- 
tation increases, as well as the effect of advection (see 
Appendix |c|). Thus, the toroidal magnetic field at the 
disc surface decreases, decreasing accordingly the mag- 
netic torque. This has two consequences. First, the mass 
load is sfightly lowered (see Eq.(^)), providing a higher 
Alfven speed at the Alfven surface and a lower fastness pa- 
rameter. Second, the jet is less powerful {(Jsm decreases) 
and the magnetic force (Fj_, see Eq.(15)) is less efficient in 
opening the jet: va/to decreases while za/to increases. In 
order to accelerate plasma up to the Alfven surface, the 
magnetic structure has to provide more power (through 

and F||), thereby reducing the available current flow- 
ing inside a given magnetic surface (gA decreases) . There- 
fore, powerful jets with current still available at the Alfven 
surface require a high diffusivity level («,„ ^1). 

We can now generalize our understanding of the 
asymptotic behaviour of super-Alfvenic jets by using loa 
instead of ^. Indeed, although the jet parameters are 
mostly described by ^, the jet asymptotic behaviour is 
drastically sensitive to u>a- 



Fig. 5. Influence of the turbulence parameter am on other 
MAES parameters (see Fig. for e = 10~^ and various ejec- 
tion indices: ^ = 0.004 (solid line), 0.005 (dotted line), 0.01 
(short-dashed line) and 0.02 (long-dashed line). The minimum 
level of MHD turbulence is limited by the value of the induced 
toroidal field that allows trans- Alfvenic jets {uja > 1). The 
maximum level is arbitrarily fixed to unity. Note that for fixed 
Um, the fastness parameter uja grows with decreasing ^ (see 
text). 



4.2.3. Influence of the magnetic diffusivity 

In the above two sections, we restricted ourselves to discs 
with am = 1- Figure ^ shows that it has a profound 
effect on the fastness parameter uja- at constant ^, the 
smaller am the smaller uja- Since trans- Alfvenic jets re- 
quire UJA > 1, we obtain that there is a minimum turbu- 
lence level required. Because usual dimensional arguments 
restrict a,„ to unity, we conclude that steady state MAES 
require an MHD turbulence with 0.1 < am 1. 



5. Asymptotic behaviour of self-similar jets 

5.1. The fate of self- similar, non-relativistic jets 

When integrated up to "inflnity", all our solutions dis- 
play the same behaviour: after an initial widening, the 
magnetic surfaces reach a maximum radius rt and then 
start to bend towards the jet axis (see Figs. ^ and 
The acceleration efficiency is very high, since the magnetic 
structure converts almost all the MHD Poynting flux into 
kinetic power (cr <C 1), allowing therefore matter to reach 
its maximum velocity (see Figs. ^ and ||). They finally stop 
at a finite distance, with a cylindrical radius larger than 
the Alfven radius va (Figs. ^ and |^). We are then faced 
to the following obvious questions: Why do our solutions 
stop ? and why do they always recoUimate ? 

The refocusing of the magnetic surfaces towards the jet 
axis stops because the flow meets the fast-magnetosonic 
(FM) critical point (see Fig. [l0| ). The associated Mach 
number is deflned as t'^ = V^/^fm^ where V ~ Ur a.t those 
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Opening angle 



Pitch angle 




Fig. 6. Poloidal magnetic field lines for e = 10 ^ and = 1, 
for ^ = 0.05, 0.04, 0.03, 0.02, 0.012, 0.01, 0.009, 0.007 and 0.005 
(the maximum radius increases with decreasing ejection index) . 
The thick line connects the position of the Alfven point for each 
solution. Note that in the range of allowed ejection indices for 
Qm — 1, the appearance of jets remains quite variable: for 
an anchoring radius ro ~ 10~^ AU from a young star, jets 
can propagate from 1 to 10^ AU of the central source, with 
a maximum radius ranging from 1 to approximately 300 AU. 
Note the logarithmic scales: small ^ jets recoUimate with angles 
smaller than one degree (Fig.^ . 

altitudes (FP95). In the limit of super- Alfvenic speeds 
((? ~ 1), this Mach number can be written as 



UJA 



/2A- 3 



1/2 ^2 



(52) 



where 9 is the jet opening angle (Fig. 0). In the other 
extreme limit {g ~ 0), one has ~ sin^ 6*. Thus, as 
the jet refocuses towards the axis (but with r > ta), its 
opening angle grows and = 1 is unavoidable. 

The integration scheme stopped because no regularity 
condition was imposed. A priori, such a critical point could 
(and should) be smoothly crossed if, inside our parameter 
range, we could obtain two characteristic behaviours at its 
vicinity. But we never obtained "breeze" -like solutions, all 
of them ending like in Fig. 10, in a way very similar to 



Fendt et al. ( |1995| ). Since we believe that our parameter 



space is weakly affected by the self-similar ansatz used, it 
seems to us doubtful that any physical trans-FM solution 
could be found. One should bear in mind that, to this date, 
there is no semi-analytical jet solution propagating to in- 
finity that crosses all three critical points (slow, Alfven 
and fast). Providing a precise answer to whether such a 
solution could exist is beyond the scope of the present 
paper. However, we will give next section an argument 
against such a possibility, based on the physics that leads 
to this critical situation. 





Log x(a) 



Fig. 7. Jet opening angle {6 — aTcsm{B,./Bz), left pannel) 
and pitch angle (arctan(— B^/Bp), right pannel) in degrees, for 
e = 10"^, a™ = 1 and ^ = 0.005 (solid line), 0.01 (dotted line), 
0.02 (short-dashed line) and 0.05 (long-dashed line). The jet 
opening angle is positive when the jet widens, negative when it 
undergoes a recoUimation. While jets with high ejection indices 
refocuse with an angle of order —10°, the others are almost 
cylindrical, with an angle smaller than 1° (but negative). 

Note that all jets displayed here achieve a super 
fast-magnetosonic poloidal speed, namely an usual fast- 
magnetosonic Mach number — u^/Ypj^j > 1. The last 
critical point (t^ = 1) would in principle constrain another 
parameter of the MAES (e.g. the ejection index ^), thus 
leaving free (but severily bracketed) mostly one parame- 
ter, the disc aspect ratio e {am is supposed to be provided 
by calculations of MHD turbulence inside the disc). 

RecoUimation of magnetic surfaces is an old result: 
all BP82's solutions displayed such a "turning radius", 
while PP92 found that such a behaviour could be generic 
for cold jets (they did not us e any self-similar ansatz). 
Contopoulos & Lovelace ( 1994 ) however, using also a self- 
similar ansatz, obtained jets with different asymptotic be- 
haviours: recoUimating, ever-widening and oscillating. But 
these solutions were obtained by varying the magnetic 
configuration parameter (/3) independently from the oth- 
ers (A and n) and so, cannot describe jets from Keplerian 
discs. These different behaviours reflect the richness of the 
equations governing MHD jets. 



Following PP92, let us define x — m rj^/r as a mea- 
sure of the jet widening. At the turning point (B* = 0), 
Bernoulli equation writes 



Xt 



(53) 



which, in the limit ~ 1, provides the cubic 

x't-{2X-3)K^X\xt~l)-0 . (54) 

In that limit, one has x = "-V2 - 1 (PP92). Solving the 
cubic and using this last expression, allows us to show that 
jets will recollimate when they reach 



> = 2kA(2A- 3)^/^ - 3 . 

In the other limiting case {gt - 
rectly the condition 



0), Eq.| 



>n1 = kA(A - 3) 



1/2 



(55) 

provides di- 
(56) 
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Fig. 8. Ratio a of the Poynting flux to the kinetic energy 
flux, logarithm of the Alfvenic Mach number m^, total cur- 
rent I flowing within the magnetic surface (normalized to the 
current Ism provided at the jet basis), and logarithm of the 
jet density along any magnetic surface, for various ejection in- 
dices. These curves were obtained for the same parameters as 
in Fig. ^. The jet density is normalized to the disc midplane 
density po = Ma{ro) /?>'ivqp,£^Q,or1, at the radius ro- The jet 
acceleration depends on how much a magnetic surface widens 
{m? = Pa/ p), which is possible only if the current decreases. 
Here, this acceleration is very efficient, the magnetic structure 
feeding plasma with almost all its power (cr <C 1). 
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Log x(a) Log x(a) 

Fig. 9. Components of the jet poloidal velocity Up and log- 
arithm of the ratio of the poloidal to the azimuthal velocity, 
measured along a magnetic surface for various ejection indices 
(see Fig. For these typical solutions, the jet always reaches 
its maximum velocity (see text), mainly as a vertical compo- 
nent. For small ejection indices, a full relativistic treatment 
should be used for jets around compact objects. Inside the 
disc, matter is being accreted with a velocity of order e the 
Keplerian velocity. The turning point [ur — 0) occurs roughly 
at the disc surface {x = 1, see FP95). 



Both criteria, similar to those found by BP82 and PP92, 
are indeed verified numerically for extreme cases (all so- 
lutions have > 1). However, we used only Bernoulli 
equation to derive them: it is the jet transverse equilib- 
rium that will tell us whether or not such a possibility is 
indeed verified. Thus, we haven't provided yet a satisfac- 
tory answer to why do jets recoUimate. 

5.2. Why do self-similar jets from Keplerian discs always 
recoUimate ? 

In order to understand what forces the jet to bend towards 
the axis, one can fruitfully look at another form of the jet 
transverse equilibrium equation (20), namely 

IJ-or 

where = Va • V/| Va| provides the gradient of a quan- 
tity perpendicular to a magnetic surface (V ±Q < for a 
quantity Q decreasing with increasing magnetic flux) and 
TZ, defined by 
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Fig. 10. Characteristic velocities in units of the disc sound 
speed Q,oho for ^ = 0.009, when e — 10~^, am ~ 1. The jet 
becomes super-SM slightly above the disc, meets the Alfven 
surface at za — 9.6ro, where ro is the field line anchoring radius 
on the disc, and finally stops when it meets the FM critical 
surface. See FP95 for the expressions of these critical velocities; 
in particular, V = n ■ u is not the jet poloidal velocity. 
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is the local curvature radius of a particular magnetic sur- 



face (Appl & Camenzind 1993a). When TZ > 0, the sur- 
face is bent outwardly while for 7?, < 0, it bends inwardly. 
The first term in Eq.(57) describes the reaction to the 
other forces of both magnetic tension due to the magnetic 
surface (with the sign of the curvature radius) and iner- 
tia of matter flowing along it (hence with opposite sign). 
The other forces are the total pressure gradient, gravity 
(which acts to close the surfaces and deccelerate the flow, 
but whose effect is already negligible at the Alfven sur- 
face), and the centrifugal outward effect competing with 
the inwards hoop-stress due to the toroidal field (Fig. [TTI ). 



f= .005 



f= .05 




Log x(a) 



Fig. 11. Forces perpendicular to a given magnetic surface, 
along it, in units of Bp/fior, for am = 1, e = 0.01 and ^ = 0.005 
(left), ^ = 0.05 (right). The solid line is the sum of gravita- 
tional, centrifugal and hoop-stress, the dotted line is the total 
pressure gradient and the dashed line is the magnetic structure 
response (namely, the term (1 — iv?)r /TV). Confining forces are 
negative. In the high-WA limit (left), it is the hoop-stress that is 
responsible for recoUimation (Log xa = 1.92, Log xt — 3.6). In 
the \ow-uja limit (right), it is the pressure gradient associated 
to the poloidal field (Log xa ~ 2.1, Log Xt — 2.2). 



For current-carrying jets (high uja, 9a ^ 1), recoUi- 
mation occurs because of the constriction effect of the 
toroidal magnetic field (as first discovered by BP82). For 
such highly super- Alfvenic solutions, the jet transverse 
equilibrium depends mostly on the balance between this 
force and the centrifugal one. Therefore, it will undergo 
recoUimation only if the jet widens enough and reaches 
the cylindrical radius 



where 



V^ln/ 
V I Inr 



/3-1 



d\nx 



(59) 



(60) 



describes how strongly the current varies across the mag- 
netic surfaces at (zt ,rt). In the expression of the turning 
radius, the other term scales as implying that for 

smaller ejection indices, the jet must open more in order 
to reach r^. However, Fig. |^ shows that rt increases with 



^ much quicker than ^ ■^/'*. Thus, Ct plays an important 
role and decreases very rapidly with ^: we found numeri- 
cally that for e = 0.01, a„ = 1, we get Ct ^ 1, 10"\ 10"'' 
for ^ — 0.015, 0.01, 0.005 respectively. In fact, Ct is more 
generally a function of the fastness parameter lua (which 
increases for decreasing ^): the faster the "rotator" and 
the bigger the maximum radius reached (see Fig. |l2|). The 
centrifugal picture is therefore quite appealing here, but 
one should not forget that uja is a measure of I a ■ 

For current-free jets (low uja, Qa <C 1), rt ^ ta and 
the flow is only slightly super- Alfvenic. Here, the ratio of 
the toroidal magnetic force to the centrifugal one writes 



V I In/ 



In 5 



Ct9t 



(61) 



with of order unity but <C 1. Thus, the force respon- 
sible for recoUimation can only be the magnetic pressure 
gradient associated with the poloidal field. Indeed, the jet 
transverse equilibrium (57) implies that this force is nega- 
tive, hence trying to close the jet. If both centrifugal force 
and magnetic force associated to the toroidal field cannot 
overcome it, the jet cannot avoid a recoUimation (Fig. pT| ). 
If UJA > 1 is a necessary condition to obtain trans- Alfvenic 
jets, it is obviously not sufficient: the "rotator" must be 
fast enough to allow the propagation of jets much farther 
away from the Alfven surface. Since this result appears 
at the Alfven surface vicinity, we believe that it is real, 
namely independent of our self-similar modelling. 




Fig. 12. Ratio of the maximum radius achieved rt to the 
Alfven radius, as a function of the fastness parameter ua, 
for am = 1 and e = 10~^ (solid line), lO"'^ (dotted line), 
10~^ (short-dashed line), 7 10~* (long-dashed line). Regular- 
ity conditions at SM and Alfven points fix the range of allowed 
UJA- The jet is extremely sensitive to this parameter: jets with 
1 < UJA < 1.5 reach only a few times the Alfven radius before 
undergoing recoUimation. 



16 Jonathan Ferreira: Magnetically-driven jets from Keplerian accretion discs 




Fig. 13. Isocontours of poloidal current density (Jp, dotted lines) and poloidal magnetic field lines {Bp, solid lines) for a cold 
jet with 5 — 0.01, launched from a disc of e = 10~^, Om = 1. The current circuit displays a butterfly-like shape, characteristic 
of tenuous ejections < 1/2). The current flows down the jet axis, enters the disc at its inner edge Vi and returns along the jet 
itself. Force-free solutions would have Jp \\ Bp, which is not quite achieved here, even for smaller values of ^. The jet carries its 
own current, building up a global electric circuit as it propagates through the interstellar medium. 



We can now try to understand the fate of cold jets 
from Keplerian discs, once they start to recoUimate. At the 
turning radius, matter has almost reached its maximum 
poloidal velocity, with 3> 1 (for high lua jets). There is 
therefore no way to slow down the poloidal motion re- 
mains too small) and, despite the decrease in the jet radius 
(leading to a decrease in rri^ = pa/ p), the poloidal veloc- 
ity remains roughly constant. Matter carries away the field 
with it, strongly decreasing the jet pitch angle. Thus, the 
total current / (or equivalently Bip) goes to zero. Eventu- 
ally, the dominant magnetic pressure effect in Eq.(57) be- 
comes the one due to the poloidal field. Because of recoUi- 
mation, the pressure gradient associated with the poloidal 
field changes its sign and pinches the jet too (some sort of 
depression). As a result, the jet is forced to bend towards 
the axis, with a curvature radius TZ becoming infinitely 
negative. This behaviour was also obtained by Contopou- 
los & Lovelace (1994) and lead these authors to propose 
that the assumption of steady-state should break down at 
these distances. More generally, any steady-state jet with 
zero poloidal current would be asymptotically parabolic 
(Heyvaerts & Norman 1989 ). Such a situation is very far 
from the actual rccoUimating regime. Thus, it seems rea- 
sonable to guess that stationarity breaks down at the jet 



"end" (where = 1 is met but not crossed) and that a 
shock is formed there. 

It is interesting to note that, althought they used 
the same self-similarity ansatz as we did, Contopoulos & 
Lovelace ( 1994| ) found also non-recoUimating (for /3 < 1) 
and oscillating (for /3 > 1) solutions. But, as already 
said, those were obtained with parameters that do not 
correspond to disc-driven jets and achieved regimes with 
n? < 1, thereby allowing a different asymptotic state. 
On the same line of thought, any neglect of the poloidal 
field (like in Contopoulos 1995) would miss the magnetic 
depression effect due to recollimation, hence modifying 
also the jet asymptotic equilibrium (and possibly allow- 
ing trans-FM solutions by fine-tunning a parameter). 



5.3. Why do self- similar jets widen so much ? 

For current-carrying jets (Fig. ^3|), it has been shown 
that jets recoUimate through the constriction action of 
the magnetic field. This is possible because the jet radius 
keeps on opening, allowing a huge acceleration efficiency 
as well as a decrease of the matter rotation rate. Why do 
we obtain such a systematic behaviour ? 
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To address this question, it is worthwhile to look at 
the Grad-Shafranov equation (20), written in the following 
form 



(1 - m^) J0 = Ja + Jk 



(62) 



It expresses that the toroidal current is generated by two 
main sources, namely 



Jx=pr{ — + {1- g)ny— + gVL,——^ 
y da da da 



(63) 



that depends mostly on how the magnetic lever arm be- 
haves from one magnetic surface to another, and 



9 Va ^, 
= m • V m p 

Mo?" 



2fio 



da 



(64) 



strongly dependent on the mass load at each surface. 
It is clear that the behaviour of global quantities has a 
tremendous importance on the jet equilibriimi, thus cast- 
ing doubts on approaches that are not truUy 2-D. In 
the self-similar description of disc-driven jets and within 
our parameter range, one always achieves a turning point 
(S* = 0) with (T( <C 1. Since the toroidal current at this 
point. 



1 dB, 
Mo dr 



« /3(2 - /?) 



(65) 



is positive, recoUimation is allowed only if the currents J* 
or are negative. At this turning point, one source writes 
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for all solutions. On the contrary, the other current source 
strongly depends on g (thus on uja)- for 3t — 1 it writes 



7* 



pnlr 



2Bo 
and for gt — 

pnlr 



o , dhiA , 
3 + 2A — 1 



din To 
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(67) 



(68) 



Self-similar solutions necessarily imply that ratios like 
rA/fo and pa/ Po are constant through all the jet. Their 
logarithmic derivatives are therefore zero and both cur- 
rent sources are indeed negative at the turning point 
(dlnpo/dlnro = 2/3 - 3 = -3/2 + £, < Q, FP93a). Thus, 
our systematic behaviour arises from both our MAES pa- 
rameter space and the fulfilment of this condition. On the 
contrary, any jet with gt ^ ^ (high uja) but satisfying 



d In r^/rp 1 
dlnro 2 
din p a/ Po ^ dlnpo 



din To 



din To 



(69) 



would never allow > and therefore, would not un- 
dergo recoUimation. This sufficient condition expresses 
that high-WA jets would not widen so much as they do 
here, the magnetic structure keeping stored a significant 
amount of the power provided by the disc. Since our pa- 
rameter space is only slightly dependent on our self-similar 
modelling, we claim that it describes realistic boundary 
conditions for any given magnetic surface. Therefore, we 
expect that any model of a magnetically-driven jet would 
also undergo recoUimation if conditions (69) are not met. 
"Over-widening" and then recoUimation would be general 
features, independent on the analytical model used, of cold 
jets from discs described by constant parameters. This is 
consistent with PP92, who did not use any self-similar as- 
sumption, but found recoUiniating solutions by making the 
approximation uja (their A parameter) constant through 
the jet. Sakurai ( |1987| ) performed full numerical simula- 
tions of jets from accretion discs that crossed the three 
usual critical points (measured with the poloidal veloc- 
ity), but without taking self-consistently into account the 
disc. Because of the limited computational domain, it is 
difficult to figure out whether or not his solutions display 
a recoUimation (see the vertical scale involved for small 
^ in Fig|6|). However, his initial split-monopole geometry 
for the magnetic configuration forbids the development of 
a jet with a constant ^ (see his Fig. 3). To determine if 
the final configuration is realistic would require a full 2-D 
treatment of the disc-jet interrelations. 

For current-free or 1ow-wa jets {gt 0), non-recoUi- 
mating solutions could in principle be obtained if 



din pa/ Po ^ din Po ^^^^ 



din Tn 



d In r, 



A 



(70) 



is verified. This necessary condition is much more strin- 
gent that in the case of high-cj^ jets, and we have strong 
doubts that it could ever be met. As already mentionned 
however, we expect that the behaviour of these kind of 
jets is mostly independent of our modelling. 

It is noteworthy that differential rotation is of so great 
importance for both jet formation and coUimation. Indeed, 
it has already been found that without the counter cur- 
rent due to the disc differential rotation, no magnetically- 
driven jet could be possible. Here, we find that in the 
absence of differential rotation (dri,/da — 0), one of the 
two sources of toroidal current remains positive for any gt, 
namely 



2Bo 



3 + gtA 4 + 2 



dlnA 
din To 



(71) 



Thus, differential rotation of magnetic surfaces has a 
tremendous importance in recoUimating outflows, or more 
generally, in their asymptotic behaviour. This is a very im- 
portant effect and shows that jets from a rigidly rotating 
object (e.g., a star) will certainly have a different asymp- 
totic behaviour than jets from accretion discs. In partic- 
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ular, they can display non-recoUimating or oscillating be- 



haviours (Tsinganos & Trussoni 1991, Sauty & Tsinganos 



1994) 



MAES are expected to be settled in the innermost part 
of a larger accretion disc, where both the energy reservoir 
and magnetic field strength are bigger. Such a picture im- 
plies that a viscous-like transport of angular momentum 
acts in the outer disc until magnetic braking due to the jet 
overcomes it in the inner regions. One would then obtain 
a transition, as Vo decreases, from ^ = to ^ ^ 0.01 with 
the corresponding changes in taIto (from infinity to a fi- 
nite value) and pa/ Po (from zero to a finite value). Thus, 
the conditions (69) could be naturally fulfilled in a self- 
consistent picture taking into account realistic boundary 
conditions. Needless to say that only 2-D numerical cal- 
culations could achieve it. 

6. Summary and conclusion 

In this paper, we have covered all dynamical processes at 
work in a stationary, weakly dissipative MAES of bipo- 
lar topology (see Appendix ^ for quadrupolar), by con- 
structing continuous solutions from the accretion disc to 
super- Alfvenic jets. We summarize below our findings and 
discuss their implications. 

(1) We have demonstrated that the usual parameters 
for magnetically-driven (cold) jets, namely the magnetic 
configuration, lever arm and mass load, are intrinsically 
related to the disc ejection efficiency ^. Thus, realistic jets 
would be described by a distribution of ^ with the mag- 
netic surfaces (or their anchoring radius). 

(2) We showed that this parameter lies in a very nar- 
row range for cold jets, namely 0.004 ^ ^ ^ 0.08. A min- 
imun ejection efiiciency is required for the disc to find a 
quasi-MHS equilibrium. A crude modelling of this sensi- 
tive equilibrium leads most probably to unstable regimes, 
the disc being too much pinched by the Lorentz force. On 
the other hand, a maximum mass load arises from the con- 
straint of accelerating matter up to super- Alfvenic speeds. 
The threshold for the minimum ejection index is raised for 
decreasing disc aspect ratio e and turbulence parameter 
am- Both constraints impose that these parameters must 
verify 0.0005 < e < 0.1 and 0.1 < a™ < 1. 

(3) As a result, the current must enter the disc at 
its inner edge, flow back in the jet and close by flow- 
ing down along the axis. This has two important conse- 
quences. First, altho ugh t he disc can afford larger ejection 
efficiencies (Ferreira 1994 ), they would be non-steady, the 
flow not being able to reach the Alfven surface. Second, 
this implies a strong influence on what's going on at the 
axis, thus suggesting a coupling with the inner central ob- 
ject magnetosphere. This work is under progress. 

(4) The existence of a minimum ejection index has 
strong implications on jet energetics. Indeed, it forbids to 
construct jet models where the mass load is arbitrarily 
small (and thus, an asymptotic velocity arbitrarily high). 



Nevertheless, in the range allowed for ^, the ratio of the 
total ejection rate (in one jet) to the accretion rate, / = 
MjejMae-, would vary between 10"'^ and 10~^, depending 
on both ^ and the radial extension re/r^ of the magnetized 
disc (Ferreira 1996| ). This is in complete agreement with 
recent estimates made by Hartigan et al. (1995), who flnd 
a typical value of 10~^ for jets from young stars. 



The j et m agnetization parameter (Michel 19691 , Ca- 
mcnzind 1987), 



shows that relativistic jets could be possible around com- 
pact objects. If Compton drag can be avoided (Phinney 



1987 ), the mean bulk Lorentz factor 7 that such jets could 



achieve, 



7-1 
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24/^ 



(73) 



strongly depends on how much power the magnetic struc- 
ture keeps stored in its asymptotic regime (see FP95 for 
the definition of rjnb). As an example, for aoo — 1 (final 
kinetic energy flux comparable to the Foynting flux, Li et 
1992), one gets 7 between 1 (^ 
= 2r. 



al 



8 (^ = 0.004, 



0.07, re = lOOr,) and 
with typical values lying between 



1 and 4. Thus moderate mean bulk Lorentz factors are 
likely to be achieved by cold jets from compact objects. 
In the context of YSOs, the mean jet velocity writes 



= ar,(l + aoo)-^/^(|^ 



1/2 



(74) 



where ^tri is the angular velocity at the inner disc radius. 
For CToo — 1, jets can reach a velocity between 1 (^ = 0.07, 
re = lOOrj) and 10 (^ = 0.004, re = 2ri) times this veloc- 
ity, with typical factors between 2 and 6. For structures 
settled at ~ SOrQ (10 stellar radii for a typical T-Tauri 
star), this provides jets with a mean velocity between 100 
and 500 kms~^, in agreement with observations. 

(5) All the above results are general to magnetically- 
driven jets from Keplerian accretion discs. We constructed 
global solutions, from the accretion disc to a super- 
Alfvenic jet, using a self-similar ansatz. It has been shown 
that the jet asymptotic behaviour strongly depends on 
the fastness parameter uja, which describes how fast is 
the magnetic "rotator" and is a measure of the amount of 
current still available at the Alfven surface. This parame- 
ter must be bigger than (but of the order of) unity and is 
very sensitive to the physical conditions inside the disc. All 
our solutions display the same behaviour: the jet widens 
until a very strong acceleration efficiency is achieved, all 
the available power being eventually transferred into ki- 
netic power; the centrifugal force decreases so much that 
the Lorentz force pinches the jet, making it to recollimate. 
The bigger uja and the larger the maximum radius reached 
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by the jet. Inside our parameter range (which is weakly 
affected by our self-similar assumption), this behaviour 
stems from having a constant ejection index ^. Therefore, 
we expect that non self-similar jets from Keplerian discs 
described with constant parameters would also display 
recoUimation, unless realistic 2-D boundary conditions are 
taken into account. Such boundaries should concern the 
transition from the outer viscous disc to the inner magne- 
tized one, as well as a possible interaction with the central 
object. 

(6) If magnetized discs are driving jets over a wide 
range of radii, then they can probably be described by 
an almost constant ejection index ^ (and self-similar solu- 
tions are not too bad an approximation). In those circum- 
stances, one would expect that these jets undergo a rec- 
oUimation and then a shock, making the whole structure 
unsteady. Whether or not such a shock is terminal or a 
"magnetic focal po int" (" Gomez de Castro & Pudritz 1993| , 
Ouyed & Pudritz 1993 ) remains to be carefully worked 
out. Besides the shock signature, such a structure could 
be detected through its apparent lack of radiation coming 
from the magnetized disc itself (FP95). However, one has 
to bear in mind that the fate of jets could be strongly 
modified by the external medium (Appl & Camenzind 
1993a). Indeed, an external confinement could forbid a 



natural "over-widening" , thus enforcing a redistribution of 
the current density inside the jet. By this way, asymptotic 
solutions with V||/ = but with non- vanishing current 
/ (in contrast with what is obtained here) could perhaps 
be achieved. Such an hypothesis should deserve further 
investigation. 

(7) Finally, a major challenge remains the question of 
the source for the required magnetic diffusivity. Indeed, 
Heyvaerts et al. ( |1996D showed that in a disc braked by 
viscous stresses, any instability (possibly magnetic), trig- 
gering a turbulence with an injection scale of order the 
disc thickness, would provide TZm ^ 1- Such a situation is 
therefore incompatible with magnetically-driven jets. As 
a consequence, these powerful jets require physical condi- 
tions that were not yet investigated in discs like those, for 
example, met at the interaction with the central object 
magnetosphere. On the other hand, thermally-driven jets 
(but magnetically confined) should be viewed as a pos- 
sible alternative. The jet power could be reduced down 
to a level comparable to the disc luminosity (see Eq.(14), 
with A ~ 1). Such a picture is appealing because it allows 
jets along with radiating discs and offers a smooth transi- 
tion between "viscous-like" discs and MAES (Ferreira, in 
preparation). 
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A. Appendix: Quadrupolar topology and the pro- 
duction of jets 

In this appendix, we briefly investigate the structure of a 
Keplerian accretion disc, thread by a large scale magnetic 
field of quadrupolar topology. Such a topology is described 
by an even toroidal field and a magnetic flux a(r, z) (see 
Eq.(||)), being an odd function of z. The magnetic field 
inside the disc follows the plasma, with B^.o = 0, Br,o < 
and Bf^ o > (subscripts "o" refer to quantities at the 
disc midplane). Since both plasma angular velocity and 
poloidal magnetic field increase towards the center (we 
don't take into account a possible boundary layer between 
the disc and the central object), the disc c an be divided 
into two distin ct reg ions (see Lovelace et al. 1987 , Khanna 
& Camenzind 1992|) . The inner region is characterized by 
a positive vertical current density at the disc midplane 
{Jz,o > 0), but negative at its surface (J+ < 0). The outer 
region displays the opposite behaviour, with J^^o < and 
J+ > (the transition between the two is of course at the 
radius where Jz^o = 0). 

One would like to build up a Keplerian disc (that is, in 
quasi-MHS equilibrium in both radial and vertical direc- 
tions), accreting towards the center {ur,o < 0) and giving 
rise to a magnetically propelled jet. This last demand re- 
quires the following (necessary) conditions: (a) open field 
lines {B^ and both positive), (b) a positive MHD 
Poynting flux (for 2: > 0, namely B'^ < 0) and (c) a MHD 
acceleration (Vy/ > 0, see Eq.(15)). 

The inner region is unable to steadily produce MHD 
jets. This arises from the topology itself, which provides 

< (therefore < 0): the magnetic structure would 
slow down any plasma that could have been (thermally) 
ejected out from the underlying disc. 

On the contrary, the outer region can in principle pro- 
duce jets. Let us then examine its structure. Since the 
Lorentz force is accelerating the matter at the disc mid- 
plane, accretion is achieved only if the viscous torque is 
dominant (A < 1). Thus, the accretion velocity is 



2(i_A)l^ 
r 



(Al) 



where i'v,o is the anomalous viscosity (measured at the 
disc midplane). Both radial and vertical components of 
the Lorentz force are positive, thus counter-acting gravity 
inside the disc. If the disc is in quasi-MHS vertical equi- 
librium, then the deviation from Keplerian rotation law is 
of the order e^. At the disc surface, the vertical Lorentz 
force changes its sign, keeping nevertheless > if 
> J^B^ /B^ is satisfied (condition (c)), with J+ > 0. 
Condition (a) is fulfilled at the disc surface only if the 
vertical scale of variation of the magnetic fiux is of the 
order of the disc scale height. Steady-state diffusion of the 
poloidal field (Eq.(4)) requires then a "poloidal" magnetic 
diffusivity such that h'm,o — suoh. The toroidal field at the 
disc surface can be written as = ^ + BY^, where 
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BJ^ is provided by the disc differential rotation (Eq.(5)). 
Since B'f must be negative (condition (b)), 



-Bro 2 Uo v'„ 



(A2) 



where v'^ ^ is the "toroidal" diffusivity (see Sect. 2.1). 
With these estimates, one obtains that with the follow- 
ing restrictions, namely 



Bl 



< 1 



(A3) 



2(1- A)£2 



where A < 1 and the plasma pressure Pq — po^gh^ (quasi- 
MHS vertical equilibrium), quadrupolar topologies fulfill 
the disc requirements as well as the jet conditions (a) and 
(b) (condition (c) requires a full calculation of the vertical 
structure). Note that above the disc surface, the current 
system is similar to the one of a bipolar magnetic topology: 
the magnetic energy that feeds the jet arises from the disc 
plasma itself, extracted where < (in a layer around 
Jz =0). Since the torque there cannot be as strong as in 
the bipolar case, one gets the general result that jets from 
quadrupolar topologies would be weaker. 

However, this configuration demands very peculiar 
conditions on the disc turbulence. Indeed, a situation giv- 
ing rise to A < 1 requires I'm {^ — -^)^'^'^m ^ith v'^ ^ i^y, 
whereas A '-^ e is achieved for ''^ £v' ~ e^v^, and 



A 



for j/„ 



Such a situation is against our 
current understanding of turbulence, where all transport 
coe fficien ts should achieve a comparable level (Pouquet et 
al. 1976 ). Therefore, it is dubious that such a topology 



could produce MHD jets. It could nevertheless play a role 
in a quadrupolar "cored apple" circulation around proto- 
stellar sources (Henriksen & Valls-Gabaud 1994, Fiege & 
Henriksen 1996[ ). 



B. Appendix: Toroidal field at the disc surface 

We start with Eq.(^o|), by making a second order Taylor 
expansion of all the quantities and neglecting the advec- 
tion term, and we obtain 



V'oJo 



(Bl) 



where F — ivo/a^v'^ is a measure of the degree of 
anisotropy of the MHD turbulence inside the disc. The 
disc radial equilibrium provides the angular velocity = 



r2fc(l +ti'), where the deviation from Keplerian rotation 
law comes mainly from the radial magnetic tension, 



7^, 



9 



(B2) 



with V = Pol p ~ 1 -I- v'^x^ /2. The above expression shows 
that, as the density decreases vertically, the magnetic ef- 
fect increases and the plasma rotates with a lower rate. 
This gives rise to an enhanced accretion velocity above 
the disc midplane (and a corresponding source of toroidal 
current, see Figs. 4, 5 and 8 in FP95). The disc vertical 
equilibrium provides the density profile, namely 



1 



pq 



2 ^2 



(B3) 



where the first term is the tidal force, the second one is 
the magnetic pressure due to shear {B^) and the last one 
is the magnetic pressure due to curvature {B^)- These 
three effects are comparable (see Fig. ||). Using that 
the "toroidal" magnetic resistivity decreases vertically as 

77m ='7o(l- 6^02)2/2), WC get 



7 

Jo 2 



with 5 = V [1 + p^] - do (B4) 



and therefore, Ba 



Bi 



B 



vn 



where Bv' ~ —qBgX 



and -BJ^ = qBo5x^/6. The toroidal field inside the disc is 
then mainly measured by B^ = —qBof{x), with f{x) = 
x(l — (5x^/6). Since the magnetic diffusivity decreases on 
a disc scale height {do of order unity), the only way to 
allow a decrease of Jr on such a scale (5 > and of or- 
der unity) is to require F ~ 1. This is the key factor for 



ejection (Ferreira & Pelletier |l993b| , FP95). This, in turn, 
implies /+ — /(I) ~ 0.5 to 1. Such an important result is 
independent of the diffusivity used, as long as it decreases 
on a disc scale height. 

C. Appendix: Minimum ejection index 

The disc is not in a perfect MHS equilibrium, but there 
is a slight motion towards the disc midplane. From mass 
conservation, we obtain that this tiny motion is 



euoii - l)x 



(CI) 



Thus, this velocity increases as ^ diminishes, describing 
the basic fact that the plasma pressure gradient is less ef- 
fective in sustaining the disc against both tidal and mag- 
netic compression. This then implies that w" decreases (see 
Appendix 1^), which leads to a decrease of i5, hence an in- 
crease of B d, a t the disc surface (/+ increases). Now, a 
look at Eq. (^) shows that a vertical equilibrium will not 
be achieved for /+ too high (how high being obtained only 
by numerical means). 

Equation (^0|), taking into accoimt the advection term, 
can be written as 



V'oJo 



1 -F 1 



l)^77i^7n^ 



X 



(C2) 
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This expression shows that inside the disc, where this ex- 
pansion is valid, the effect of advcction is to decrease the 
radial current density (although this contribution, of or- 
der a^e, is negligible). On the contrary, above the disc 
when Uz > 0, this term changes its sign and contributes 
to maintain > 0. Eventually, as one goes from the resis- 
tive disc to the ideal MHD jet, this term will completely 
balance the effect of differential rotation. What the above 
expression shows, is that its influence will start at lower 
altitudes for higher am and e. Thus, to the highest values 
of both am and e corresponds the highest value of /+, 
henceforth the smallest ejection index. 
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